Paperfolding morphisms, planefilling 
curves, and fractal tiles 

Michel Dekking 
January 7, 2011 



Michel Dekking, Delft Institute of Applied Mathematics, 
Technical University of Delft, The Netherlands p.M.Dokki 



iiSQ'tudclft.nl 



Abstract 

An interesting class of automatic sequences emerges from iterated 
paperfolding. The sequences generate curves in the plane with an al- 
most periodic structure. We generalize the results obtained by Davis 
and Knuth on the self-avoiding and planefilling properties of these 
curves, giving simple geometric criteria for a complete classification. 
Finally, we show how the automatic structure of the sequences leads 
to self-similarity of the curves, which turns the planefilling curves in a 
scaling limit into fractal tiles. For some of these tiles we give a partic- 
ularly simple formula for the Hausdorff dimension of their boundary. 

Keywords: Automatic sequences, iterated paperfolding, L-systems, IFS. 

1 Introduction 



Folding starts with folding in two. In the early sixties John E. Heighway 
discovered that the pattern of creases that results from repeatedly folding a 
piece of paper in half generates an intriguing curve in the plane when the 
creases are unfolded to 90° angles. This curve is known as Heighway's dragon 
curve. 
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Figure 1: Folding a strip of paper 1, 2, and 3 times, coding the folds. 



The abstract sequence of the symbols D and U, coding the 'down' folds 
and 'up' folds (cf. Figure [TJ is known as the paperfolding sequence^, and is 
extensively treated in Allouche and Shallit's monograph pp. 

On the monoid {D, U}* we have the morphism tt which permutes letters: 

tt(D) = U, ir(U) = D, 

and the anti-morphism r which reverses words 

t(wiW 2 . . -w m ) =w m .. .w 2 Wx for all w 1 w 2 ■ ■ ■ w m G {D,U}*. 

The process of iterated folding is allied to the composition of these two mor- 
phisms, the morphism 7rr. We therefore shortly write 

S := itt(S) for any word S in {D, U}*. 

We next define the folding convolution S * T of two words S = ai . . . a s _i 
and T — b\. . . b t -\ by 

c t / SbiSb 2 S . . . Sb t -\S if t is even , . 

\ Sb{Sb 2 S . . . Sb t ^S if t is odd. ^ } 

It is quickly checked that this operation is associative: 

(R* S)*T = R*(S * T) for all R,S,Te {D, U}*. 

Thus n-fold (sic) convolutions S* n are unambiguously defined, and since 
S< n+ V s t ar t s with S* n for all n we can also talk about the infinite folding 
sequence S*°°. 



1 We follow the convention of pQ: when A is an alphabet, elements from A* will be 
called words or strings, elements from A°° will be called infinite words or sequences. 



2 



Note that the three strings in Figure [T] are respectively D, D * D, and 
D * D * D. It can be seen in general that folding first according to T and 
then according to S gives the pattern coded by S * T. 

We end this section with a few historic remarks. The Highway dragon was 
brought to a larger audience by Martin Gardner in 1967 Mathematical 
proofs of its properties were given by Davis and Knuth three years later in [3] 
(reprinted in 2010 in p2]). Around 1975 the present author generalized the 
notion of folding by introducing the folding convolution in Equation (JTJ) . He 
announced this in two letters send to Martin Gardner. These were found by 
Donald Knuth last year in the Stanford archive. Knuth incited the author 
to write down the results with proofs. Some of these results have been 
announced in the report [5] and the paper [6] jointly written with Mendes 
France and van der Poorten. 

2 Folding morphisms 

Since the main interest in folding sequences is in the curves they generate in 
the (complex) plane, it is more convenient to describe these curves directly. 
We do this by passing from an alphabet of folding instructions {U,D} to 
an alphabet A = {a, b, c, d} of four letters, representing line segments in 
the positive real, positive imaginary, negative real, and negative imaginary 
direction. 

It does, indeed, appear that the process of iterated folding and unfold- 
ing to 90° angles can be described by iterating a uniform morphism on the 
alphabet {a, b, c, d}. For 2-folding this morphism 9 is given by 

9(a) = ab, 9(b) = cb, 9(c) = cd, 9(d) = ad. 

How can this be seen? 

Define the morphism a on the monoid A* which cyclically permutes the 
letters: 

a (a) = b, er(fe) = c, a(c) = d, a(d) = a. 
Also let the anti-morphism r, which reverses the words from A*, be as before. 

The crux of folding morphisms is that they commute with the morphism 
ar. For the 2-folding morphism this holds for the letters: 

9ar(a) = 9(b) = cb = ar(ab) = crr9(a), 
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similarly for the other three letters, and then for any word. 

Suppose now that 9 n (a) corresponds to the curve obtained by folding in 
two n times. Then if we unfold a strip folded n + 1 times, the first 2 n letters 
are equal to those we obtained before, and the last 2™ letters are equal to the 
first 2 n , but in reverse order and cyclically permuted (because of the 90° angle 
that we make half way). In other words, this word should be 9 n (a)ar9 n (a). 

But 9 n also commutes with err, hence 

9 n {a)aT9 n {a) = 9 n (a)9 n aT(a) = 9 n (a)9 n (b) = 9 n (ab) = 9 n+1 (a), 

so that indeed it follows by induction that the morphism 9 describes the 
iterated 2-folding. 

We define in general a folding morphism as any morphism 9 on A* that 
commutes with or. Such morphisms are of course determined by their value 
at the letter a, and we will require that 9(a) starts with the word oh. 

We then have a one-to-one correspondence between strings S — D ... of 
length m— 1 and folding m-morphisms. To see this, define the homomorphism 
w : {D,U}* ^ Z by 

w(D) = +1, w(U) = -1. (2) 

Then note that if S — a\ . . . a m _i, and 9(a) = eo . . . e m _i, then eo = a and 

e i = <r w ( ai - ai) (a), for i = l,...,m-l, (3) 

where a is the cyclic permutation morphism on {a, b, c, a 1 }*. 
Conversely, define the 2-block map p by 

p(ea(e))=D, p(a(e)e) = U for e G A 

Then p(e ■ ■ ■ e m _i) = p(e e 1 )p(e 1 e 2 ) . . . p(e m - 2 e m -i) gives the string back. 

Our one-to-one correspondence is an operation-preserving isomorphism 
if we consider the convolution operation * on strings and composition of 
morphisms: 

Lemma 1. Let S = D ... and T = D ... be two strings with associated 
folding morphisms 9 S and 9 T . Then 

9s*t = 9s9 T - 
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Proof: We give a purely geometrical argument: S*T describes the segments 
in the curve associated to first folding according to T, and then to S. But the 
same curve is described by replacing the segments in direction i k (k = 0, 1, 2 
or 3) of the curve described by 9t{o) by the (scaled) curves associated to 
(ar)%(a). 

The algebraically inclined reader may provide an alternative proof, using 
equations (J2J) and (j3J). □ 

In the sequel all strings S from {D, U}* will start with the letter D, and 
9 will always be the morphism on A* = {a, b, c, d} associated to S — and so 
9(a) will always start with ah. 

3 Representations of morphisms 

We now do a more elaborate analysis of the situation based on ([I]). 
As usual we will write for a set B C C, and r, t G C: 

rB = {rx : x G B} and B + t = {x + t:xE B}, 

and we will write segments in the complex plane as 

[Hi y'} = {x ^ C : x = ay + (1 — a)y', < a < 1}. 

We introduce a homomorphism / : A* — > C by 

/(a) = 1, /(&) = z, /(c) = -1, /(d) = -z, (4) 

where the operation on A* is concatenation, and that on C addition. To 
associate curves to words we define a map K[] : A* -> 2 C , by first putting 

tf[a] = [0,1], = [0,z], K[c] = [0,-1], K[d] = [0,-z], 

and then extending this by putting for any two words v, w G A* 

K[vw] = K[v] U (K[w\ + f(v)). (5) 

So for example K[abcd\ = K[a] U (K[b] + 1) U (K[c] + 1+z) U (lf[d] + z) is 
nothing else but the unit square in C. 
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For some morphisms 9 we will have a commuting diagram 



A" 



-> A* 



where Lg is multiplication by a complex number. It is not hard to see that 
Lg exists for any folding morphism and is multiplication by z = f(9(a)). 
For example for the 2-folding morphism determined by 9(a) = ab, Lg is 
multiplication by 1 + i. As another example: L a is multiplication by i for 
the cyclic permutation morphism a. 

Folding convolutions behave nicely with respect to the embedding in C: 

Lemma 2. Let S and T be two folding strings, whose curves end in zs and 
Zt- Then the curve of S *T ends in z$zt- 

Proof: We have that Lg s is multiplication by zs = f(9s(a)), and Lg T is 
multiplication by zt = /(^t(o))- Lemma [T] implies that the curve of S * T 
ends in the point f(9s9T{a))- The result thus follows from a look at the 
following commutative diagram, recalling that f(a) = 1. 



A* 



C 



> A* 



6s 



> A* 



C 



□ 

In Sections [6] and M we will need a map from unions of Gaussian integer 
segments to itself that describes the action of 9. We consider neigbouring 
Gaussian integers z and z', defining the A-code e Z}Z i of the segment [z, z'\ by 

62,2' = e if z' — z = f(e), where e = a,b, c, or d. 

If [z, z'\ is such a segment then we define the K-9-image of [z, z'\ by 

K[9]{[z,z'}) = K[9{e z , z ,)} + Lg(z). (6) 

For example, if 9 is the 2-folding morphism, then 

K[9]([l + i,i}) = K[8(c)] + L g (l + i) = K[cd} + 2i 

= (K[c] U [K[d] - 1)) + 2i = [2i, 2i - 1] U [2i - 1, i - 1]. 
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The map [6 1 ] ( • ) is extended to unions of k segments, in the natural way: 
K[9){[ Zl ,z' x ] U...U [z k ,z' k ]) = K[9]([z u z[]) U...U K[9]([z k , z' k }). (7) 

Note that if the union Q of the segments can be written as Q = K[w] for 
some word w from A k , then by (JSJ) simply K[6](Q) = K[8(w)}. 

4 Symmetric folds 

We call a string of IPs and D's symmetric if S — S. For example the string 
DDUU coding a 5-fold is symmetric. Obviously only strings of folds of odd 
length can be symmetric. We call a morphism 9 palindromic if it commutes 
with the reversal map, i.e., if t9 = dr. 

Lemma 3. A string is symmetric if and only if the associated folding mor- 
phism is palindromic. 

Proof: Suppose S = a\ . . . a m _i is symmetric. It is easy to see that we only 
have to check that the associated 9(a) = eo ■ ■ ■ e m -\ is a palindrome word. 
First note that the symmetry of S implies that the number of C/'s in S equals 
the number of D's, and hence that w(S) = 0, where w is the homomorphism 
from (j2j). But then we have for j = 1, . . . , m — 1 

w(cti . . . a m -j-i) = -w(a m _j . . . a m -i) = w(a l . . . 

since symmetry is equivalent to w(a m _i) = —w(a,i) for % = 1, . . . m — 1. Then 
we have by Equation ([3]) that 

e m _ j _ 1 = < 7 «'(«i"-«m-i-0( o ) = a w( ai ... aj )^ =e . for ■ = 1; m _ 1; 

which exactly means that 9(a) = e$ . . . e m _i is a palindrome. 
Conversely, if 9(a) is a palindrome, then for i — 1, . . . , m — 1 

O-i p(&i— p(c m _j6 m _ j_i ) p(e m — i— \C m — i) Q, m _j, 

and thus S is symmetric. □ 
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5 Self-avoiding folding sequences 



We call a folding string S self- avoiding if the curve associated to the sequence 
S*°° does not traverse a segment twice. It will appear that all information 
that we need to decide whether this infinite sequence satisfies this property 
is in the finite word 8(abcd), if 9 is the morphism associated to S. 

It is important to observe that all words generated by a folding morphism 
have the property that letters from {a, c} alternate with letters from {b, d}. 
We will call such words D 2 -words, for reasons that will become clear soon 
(see also Figure [5]). 

For the following development it is convenient to consider curves 
obtained from K[w] by rounding of the corners (see Figure [2]). 



Figure 2: The curves K[w] and K[w] with w = 9 e (a), where 9(a) = ab . 

The point of this is that the curves K[w] do not pass through the Gaussian 
integers u + iv, with m,ugZ. In this way, if w is a © 2 -word, and K[w] does 
not traverse a segment twice, then K[w] is a simple curve in the topological 
sense (i.e., the image of an injective map from an interval to the plane). 

The central notion in the sequel will be that of a 8-loop defined for any 
morphism 9 as the set K [9 (abed)]. See Figure [3] and Figure H] for examples. 





Figure 3: The #-loop of 9(a) = ababad. 
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Lemma 4. Let 6 be a folding morphism. Then its 6-loop is a closed curve. 

Proof: The endpoint of K[6(abcd)] is given by 

f(6(abcd)) = L e (f(abcd)) = L e ( f(a) +f(b) +f(c) +f(d) ) = 0. □ 

Theorem 1. Let S be a folding string. Then S is a self-avoiding string if 
and only if the 6-loop of S is a simple curve. 

Proof: Part 1: Suppose that the #-loop is not simple, i.e., there is a segment 
that is traced twice. We have to show that there exist n > 1 such that 
K[6 n (a)} is not simple. 

Let z = f(8(a)). The #-loop (ignoring the rounding off) consists of the 
four curves K a := K[6(a)\,K b := K[9(b)] + z,K c := K[6(c)} + (1 + i)z, and 
Kd := K[9(d)] + i z, which can be mapped to each other by a rotation and 
a translation (see Lemma |8]). Therefore, if one K e has a segment that it is 
traced by another segment of the same K e , then all four curves have such a 
segment. So the curve K[8(a)] is already not simple. 

The next possibility is that a segment of a K e is traced by K a ^, and 
again, if this happens for one e, then it happens for all e G A. Thus, since 
6(a) starts with ab = aa(a), the curve K[6 2 (a)} will not be simple. 

The last (note that K e has a segment in common with K a 3( e \ iff K a i e \ has 
a segment in common with K e ) possibility is that a segment of a K e is traced 
by i^ (J 2( e ), and now the situation is different, since this might happen for K a 
and K c , but not for Kb and Kd (cf. Figure H]), or vice versa. In the first case 
this overlap will occur in K[6(abc)}, in the second case in K[6(bad)]. So our 
goal is to show that both abc and bad will occur in some 6 n (a). 




Figure 4: The #-loop of 9(a) = abcbabadadcda. 

Now, since 6(a) starts with ab, there exists a k > 1 such that either 
6(a) = (ab) k or (ab) k a, or 6(a) starts with (ab) k c or with (ab) k ad. The 
former case can not occur, since clearly the (9-loops of (ab) k or (ab) k a are 
simple for any k > 1. In the second case 6(a) contains the word abc, or the 
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Figure 5: The directed graph B 2 . 



word bad. But since ra(bad) = abc and (to-) 3 (abc) = bad, 9 3 (a) will contain 
both abc and bad. It follows that K[8 3 (a)] will not be simple. 

Part 2: Suppose the ^-loop of S is a simple curve. The crucial observation 
here is that in whatever contorted way a curve K[8(w)] (with w a ©2-word, 
and 9 a folding morphism) lies in the plane, it will always be a subset of a very 
regular checkerboard pattern of (9-loops and anti-^-loops (ignoring rounding 
off). Here we define the anti-8-loop as the curve K[6(dcba)}. 

The regular pattern is created by the structure of the directed graph 
D 2 , which as vertices has the Gaussian integers, and directed edges between 
neighbouring vertices as in Figure [5j To avoid cumbersome notation we shall 
often identify the vertices with the corresponding Gausssian integers (which 
we will also call grid points), and the directed edges with the corresponding 
segments in the plane. 

We call a Gaussian integer z = u + iv even if u + v is even, and odd if if 
u + v is odd. 

The regions bounded by the curves K[abcd] and K[dcba] and their transla- 
tions over even Gaussian integers correspond to the white, respectively, gray 
squares in Figure [51 It follows that the curves K[9(abcd)] and K[9(dcba)} 
are #-loops (white interior) and anti-#-loops (gray interior) as in the example 
in Figure [6] (we see in fact, in this figure part of the i^-#-image i^[#](©2) 
of ©2)- Let w be a ©2-word such that K[w] is self-avoiding. Since the 9- 
loops and anti-#-loops are self-avoiding, the same will be true for the curve 
K[9{w)], which will be a subset of the full 6*-loops and anti-6*-loops pattern 
(after rounding off in the appropriate way). It thus follows by iteration of 
this observation that all K[9 n (a)] will be self-avoiding, and so by definition 
S is self-avoiding. □ 
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Figure 6: #-loops and anti-0-loops when 0(a) = ababad. 

6 Plane-filling folding sequences 

A surprising property of the curve generated by 2-folding is that it is not only 
self-avoiding, but very compactly traverses large parts of the set of Gaussian 
integers ([HI [3]). 

We call a folding string S planefilling if the curve associated to the infinite 
sequence 5**°° passes twice through every grid point lying in some arbitrary 
large circle. 

It will be convenient to reformulate this to the equivalent property that 
this curve passes through all the segments connecting Gaussian integers form- 
ing a (2k + l)-cross in D2 for every k = 1,2.... Here the (2k + l)-crosses 
are translations of (2k + l)-crosses centered at (1 + i)/2 defined as the set of 
edges in D2 lying within the square through the points k + 1 + q, (k+l)i + q, 
—k — l + q and — (k + l)i + q, where q = (l+i)/2. Note that as directed graphs 
there are actually two types of (2k + l)-crosses: those where the center loop 
runs clockwise, and those where the center loop runs anti-clockwise. In the 
sequel we will not make this distinction explicitly. 

The goal of this section is to give a characterization of self-avoiding plane- 
filling folding strings. In the next section we will see that there exist many 
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Figure 7: The 3-cross and 5-cross graphs. 



self-avoiding planefilling folding strings. 

We call a string S maximally simple, if its #-loop is simple and if the 
region the #-loop bounds (by Jordan's Theorem), only contains grid points 
that are passed twice by the 6*-loop. 

Lemma 5. Let S and T be two maximally simple folding strings. Then S*T 
is also maximally simple. 

Proof: This is similar to the proof of Lemma [TJ 

iz T (1 + i)z T 




Figure 8: The #-loops of the folds S = DUDD, T = DUU, and S * T. 



We obtain the #s*T-loop by replacing scaled unit-loops in the #T-loop by 
#s-loops. Grid points in the interior of the ^^-loop will either lie in these 
#s-loops, or will be in the interior of the scaled ^j-loop, cf. Figure |SJ In 
both cases these points are visited twice, and so the ^T-loop is maximally 
simple. □ 
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Lemma [5] forms the basis for the planefilling property. It implies that for 
a given string S with maximally simple #-loop all its # n -loops will also be 
maximally simple. However, there are some subtleties in actually proving 
that this implies that S is planefilling. 

We call a word w from A* a loop word if it is a D2-word, and f(w) = 0, 
where / is the homomorphism from Equation (j4j). 

We call a loop word maximally simple if K[w] is a self- avoiding curve that 
visits each grid point in its interior twice. Examples are the 8 words 

a k (abcd), and a k (adcb) for k = 0, 1, 2, 3, (8) 

which we call square words. 

Lemma 6. Let w be a maximally simple loop word. Then w contains at least 
one of the square words. 




Figure 9: The curve K[w] of a maximally simple loop word of length m = 40. 

Proof: A point u + iv with u = N + 1/2 and v = M + 1/2, where N and M 
are integers will be called a mid grid point. Maximally simple loop words w 
of length m are associated to particularly nice curves X[u>]. The curve K[w] 
contains exactly m/4 mid grid points in its interior. Consider the graph Q 
whose vertices are these mid grid points, and where there is an edge if two 
vertices are neighbors (i.e., the corresponding mid grid points are at distance 
\/2), see also Figure If m = 4, then w is a square word, so consider m > 8. 
Since K[w] is simple, Q is connected, and since K[w] is maximally simple, 
Q has no circuits, and so Q is a tree. A mid grid point corresponding to 
any leaf of this tree is the center of a K[a k (abcd)] or a K[o~ k (adcb)] (ignoring 
rounding of). Since m > 8 the tree has at least two leaves. At least one 
of the corresponding two square words will occur in w (one has to take into 
account the possibility that one of the two corresponding mid grid points is 
the point (1 + i)/2). □ 
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Lemma 7. Let a folding morphism 6 have a maximally simple 6 -loop. Then 
a square word will occur in 6 4 (a). 

Proof: If m = 2, then 8(a) = ab, and it is easily checked that w SQ = bcda 
occurs in 9 4 (a), so suppose m > 3. From Lemma[6]we know that w = 8(abcd) 
contains a square word w SQ . As m > 3, w SQ has to occur in at least one of 
6(ab) or 9 (be), or 6(cd). 

Obviously ab occurs in 6 3 (a) = ab. . .. We claim that also be occurs in 
6 3 (a). To see this, note that if be occurs in 6(a), then it occurs in 8 2 (a) and in 
6 3 (a), since 6(b) ends in b and 6(c) starts with c. Now either 6(a) = abc. . ., 
and then be will also appear in 6 3 (a), or 6(a) = aba. . ., and then be will 
appear in 6 2 (a) (and as argued before, in 6 3 (a)), since in that case 6(b) = 
ar6(a) ends in bcb. Finally, cd occurs in 6 3 (a), since c occurs in 6(b), and 
0(c) = (ar) 2 6(a) = cd . . .. ' 

Since ab, be and cd all occur in 6 3 (a), w SQ will occur in 6 4 (a). □ 

Theorem 2. Let S be a self- avoiding folding string. Then S is a planefilling 
string if and only if the 6-loop of S is maximally simple. 

Proof: Part 1: If the (9- loop is not maximally simple, then at least one of the 
grid points in its interior is visited only once. It follows that a (2k + l)-cross 
can nowhere occur if k is large enough. 

Part 2: Note that the 8 square words have two different types, those as- 
sociated to a ^-loop: abed, bcda, cdab, and dabc, and the other four, which are 
associated to an anti-#-loop (cf. the white, respectively gray squares in Fig- 
ure [5]). In terms of curves this means that the #-loop runs counter clockwise, 
and the anti-#-loop clockwise. 

Let w SQ be the square word occurring in 6 4 (a) according to Lemma [71 
Suppose that w SQ is of the second (clockwise) type (otherwise we replace 
below everywhere the # 4 -loop by the anti-# 4 -loop). Now consider the curve 
K[6 4 (a)} in the 6* 4 -loop. Since the type of this loop is opposite to the type of 
w SQ , the 3 grid points passed by i^[w SQ ] must be interior points of this loop, 
cf. Figure [ini left. This implies that the 10 segments indicated in Figure fTUl 
right must all be part of the (9 4 -loop. In other words: the 6* 4 -loop contains 
a 3-cross. Since w SQ occurs in 6 4 (a) this implies that K[6 8 (a)] contains a 
3-cross. 
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Figure 10: Left: K [w SQ ] in # 4 -loop. Right: forced segments in the 6* 4 -loop. 



With a figure-8 configuration we denote the subgraph of ©2 
with 6 vertices and 7 edges pictured on the right. 1 

Let Qb = [0, 1] U [1, 1 + i] U . . . U [—i, 0] be the corresponding set ~ l ^ x - 
of segments. 

Now consider K[6 8 ]([0, 1]) = K[8 s (a)] in K[P]{Q B ), see Figure EJ Since 




Figure 11: A 3-cross in a figure eight configuration at level 8. 

all the segments of i^[0 8 ]([0, 1]) lie in the interior of the (rounded off) closed 
curve corresponding to K[9 8 ](Qb) \K[8 8 ]([0, 1]), the eight boundary points of 
the 3-cross are interior points of this curve, and so it must contain a 5-cross. 
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Now note that K[9 8 ](Q B ) (rounded off) appears in K[9 w (a)), since K[9 8 (a)] 
contains a figure-8 configuration (as a subgraph of the 3-cross, both for the 
clockwise and for the anti-clockwise type). Thus K[9 16 (a)} has to contain a 
5-cross. The same type of argument yields that K[9 24 (a)} contains a 7-cross, 
and more generally that K[9 8k (a)] will contain a (2k + l)-cross, for each 
natural number k. □ 

Theorem 3. Let S and T be two self- avoiding planefilling folding strings. 
Then S *T is a self-avoiding planefilling folding string. 

Proof: This follows directly from Lemma [5] and Theorem [21 □ 

7 Existence of planefilling folds 

Let z = f(8(a)) be the endpoint of the curve associated to a m-fold. By 
Lemma [8] in Section [8]the sets K[9(a)} and K[9(c)] + (1 + i)z are congruent, 
and so are K[9(b)} + z and K[9(d)\ + iz. This implies that the interior of the 
#-loop and the square through 0, z, (1 + i)z and iz have the same area. For a 
maximally simple #-loop this area is equal to m (ignoring second order effects 
caused by rounding off), since all 4m segments of such a 9-loop contribute 
one side to a unit square. It thus follows from Theorem [2] that self-avoiding 
planefilling m- folds can only exist if m = \z\ 2 is a sum of two squares. 

Theorem 4. For any integer m > 1 which is a sum of two squares there 
exists a self-avoiding and planefilling folding string of length m. 

Proof: We will prove in fact more, namely that for every Gaussian integer 
z = u + iv with u 2 + v 2 > 1 there is a self-avoiding planefilling folding string 
whose curve ends in z. By Theorem [3] and Lemma [2j it suffices to prove 
this for those z = u + iv which are Gaussian prime numbers. Moreover, we 
can assume, by rotating if necessary, that u > 0,v > 0. We already know 
that there exists a planefilling folding string to the Gaussian prime number 
z = 1 + i. Moreover, z = u + iv is certainly not prime if both u and v are 
even, and since 



z is also not prime if u and v are both odd (and not both equal to 1). So we 
may assume in the sequel that z = u + iv is an odd Gaussian integer, i.e., 
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u + v is odd. Moreover, by exchanging u and v if necessary, we can assume 
that u is odd and v even. 

First consider the case v 7^ 0. 
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Figure 12: The octagon Q associated to z = u + iv when u = 7 and v = 4. 

We will define an octagon Jl, based on the mid grid point s = \u{l + i) 
and the grid point t = \v{l + i) (see Figure [12]). The octagon is given by its 
vertices: 

O ~ (0, s, s + it, s + (1 + -2, — «s + (1 + i)t, —is + t, —is). 

The regular way in which Q is constructed ensures that it is simply connected. 
Now consider the finite directed subgraph of D2 whose edges lie in Q. It is 
easily checked that all its vertices have in-degree 1 and out-degree 1 or in- 
degree 2 and out-degree 2, except the vertices and z, which have out-degree 
1, respectively in-degree 1 (here we use again that u + v is odd). 

Now remember that an directed Euler path is a path along a connected 
directed graph that connects all the vertices and that traverses every edge of 
the graph only once, and that the degree conditions above yield the existence 



17 



of a directed Euler path from to z ([9])- But such a path is nothing else 
than the curve given by a self-avoiding folding string S of length u 2 +v 2 , since 
the octagon Q has half the area of the square H with vertices 0, z, (1 + 
and iz. Moreover, the curve obtained by joining the four curves produced by 
rotating thrice over 90° around the center s + it = (l+i)z/2 of H is the #-loop 
of S, and it visits each of its internal points exactly twice (cf . Figure [T3|) . It 
follows from Theorem [2] that the string S generates a planefilling curve. 
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Figure 13: Four octagons forming a 0-loop region. 

When v — the octagon f2 degenerates into a square, but the arguments 
will remain the same. □ 

It can be seen that actually one can find a symmetric fold to u + iv for any 
odd Gaussian integer. The idea here is that Euler's theorem will guarantee 
a D2-word w of length (m — l)/2 with its curve filling half of Q after the 
central mid-grid square has been removed. The symmetric m-fold will then 
be given by wa t(w) if u + v — 1 is divisible by 4, or else by w ct(w). 

8 Perfect sequences: the carrousel theorem 

A string S is called perfect if it is self-avoiding, and the four curves K[S*°°], 
iK[S*°°}, -KIS* 00 ], and -iK[S*°°} together cover every edge of D 2 exactly 
once. Davis and Knuth prove in [3] that the 2-fold S = D has this property, 



18 



stating that this is "probably the most interesting property of the dragon 
curve". See Figure UM 




Figure 14: The carrousel of order 6 for S = D. 

Actually it is remarkable that a fold S can be perfect, since the set above, 
which we call the carrousel associated to S, is incompatible with the directed 
graph structure of © 2 — note, e.g., that the out-degree of the origin is four 
for the carrousel. Nevertheless, perfect folding sequences abound, and the 
following characterization holds. 

Theorem 5. Let S be a folding string that generates a self-avoiding plane- 
filling curve. Then S is perfect if and only if the curve K[9(a)] of S does not 
traverse the segment from (1 + i)z — % — 1 to (1 + i)z — %, where z = f(8(a)) 
is the endpoint of the curve of S . 

We first present a lemma which resolves the incompatibility observed 
above. 

Lemma 8. Let S be a folding string with z = f(6(a)). Then 
iK[d(a)} = K[B{d))+iz, i 2 K[d(a)} = K[6(c)], i 3 K[9(a)] = K[9(b)]+i 3 z. 

Proof: Note that these are equalities of sets, which we prove by going 
through the segments making up the sets in a certain order. To see that 
the first one holds, traverse iK[6(a)] = K[a8(a)} from iz to 0. This amounts 
to applying cr 2 r to the word a8(a) (since we reversed the direction the seg- 
ments in the curve are rotated over 180°). But since a and r commute, and 
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by definition of the folding morphism we have 



a 2 r(a9(a)) = (ar) 3 9(a) = 9(d). 
For the second equality, note that i 2 K[6(a)] = K[a 2 9(a)], and that 

a 2 9(a) = (ar) 2 9(a) = 6(c). 

The third equality follows as the first one. □ 

The segment [1 + will play a crucial role in the proof of Theorem |5j 

Lemma 9. Let S be a self-avoiding planefilling string. Then the 6-loop of S 
contains the segment [1 +i, i]. 

Proof: The statement is true when m = 2, and for m > 3 we either have 
9(a) = abc . . . , and then K[9(a)] already contains [1 + i, i], or 9(a) = aba . . . , 
and then the point 1 + % is in the interior of the 0-loop, what implies that 
the segment [1 + i, i] must be traversed by it. □ 

Proof of Theorem [5} 

Let 9 be the morphism associated to S, with 9(a) = ab. . .. First we will 
translate the mysterious condition on the curve K[9(a)] of S to a slightly less 
mysterious condition on the #-loop of S. Note that (by LemmaE]) the K[9(c)} 
part of the ^-loop (ignoring the rounding off) can be written as K[9(c)] + 
(l+i)z = -K[9(a)] + (l + i)z, where z = f(9(a)) is the endpoint of If [0(a)]. 
It follows that K[9(a)} does go through the segment 
from (l+i)z — i — 1 to (l + i)z — i if and only if K [0(c)] + 
(1 + i)z goes through the segment from 1 + % to i (see 
also the figure on the right). 

In the sequel, let 



Q n := K[9 n (a)} U iK[9 n (a)} U -K [9 n (a)\ U -iK[9 n (a)] 
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be the carrousel of order n. Since we can replace the 

2-fold D by the 4-fold D*D = DDU, we will assume henceforth that m > 3. 

Part 1: Suppose [1 + is traversed by K[9(c)} + (1 + i)z in the 0-loop. 
t^j!-K Then 9(a) must start with aba, so 9(d) must end in dad, and 
in Qi the segment [1 + i,i] is absent (cf. figure on the left). In 
1 Qi, the segment [1 + is still missing, since it could only be 
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traversed by if [0(c)] + (1 + i)z, which is missing from the incomplete #-loop 
with i as one of its internal points. Continuing in this fashion we see that 
the segment [1 + i, i] will be missing in all Q n , and so not all edges of H> 2 are 
visited. 

Part 2: We now suppose [1 + is not traversed by K[8(c)] + (1 + i)z in 
the #-loop. Note that either 8(a) starts with abc, or with aba. The latter 
case will be considered at the end of the proof. If 8(a) starts with abc, then 
the 3 rd segment of K[8(a)] traverses [1 + i, i]. Since K[8(a)] does not end in 
the point i, it will then also traverse the segment [i, 2i]. It thus follows from 
rotational symmetry that Qi will contain a 1-star, see Figure [15j 




s 



Si 



s 2 



Figure 15: The 0-star, 1-star and 2-star graphs. 



In general the k-star graph is the subgraph S k of D 2 consisting of all the 
vertices and edges lying in the square through k + 1, (k + —k — 1 and 
— (k+l)i. It will be important later that we have the following property of the 
star graphs: one obtains Sk+i from Sk by adding all edges (and corresponding 
vertices) to the vertices that are not yet of degree 4 in S^. The edges and 
vertices that are added, except for the outer 4, will be called the boundary 
B k of Sk (see also Figure [16]) . 



So 



Si 



s 2 



Figure 16: The boundaries in the 1-star and 2-star graphs. 



Now note that B\ is a closed simple loop in Si, which has all the segments 
of Sq in its interior, and that we know that Qi contains S\. It follows (recall 
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Equation ([6]) and ([7])) that Q2 contains the K-9- image of Bi, which will be a 
closed (simple) loop, that contains Si in its interior (except for the 8 extremal 
points), since the K-9- image of So is Q\. Now we get to the main idea of 
the argument: since the #-loop is maximally simple, all grid points within 
the K-9- image of B\ will be visited twice by Q 2 - But then Si has to grow^ 
till S2, and so Q2 will contain S 2 . Continuing in this fashion we find that Qk 
contains S& for all k, and so the Qk increase to cover all of Z 2 . See Figure [T71 
for an example with a self-avoiding and planefilling 8-fold. 




Figure 17: The case 6(a) = abcbadab; the carrousel Q\ and Q2 with the 1-star. 
On the right the K-6-image of B\ is indicated in orange. 

We still have to consider the possibility that 9(a) starts with aba. By 
Lemma [9] one of the appropriate translates of K[9(e)}, e = a, b or d traverses 
[1 + First, suppose that [1 + i, i] is in K[9(a)] or K[9(d)} + iz. Then 
obviously Q± contains the 1-star graph. 

Second, if [1 + i, i] is in K[9(b)} + z, then the segment [1 + i, i] will not be 
in Q\. However, since 9(a) starts with aba, 9 2 (a) starts with 9(a)9(b), and 
thus the segment + i] will be contained in Q 2 - In all cases Q2 will contain 
the 1-star graph. Thus for the stringS 1 * S the carrousel Qi will contain the 
1-star graph, and we can copy the proof of the abc case. □ 

Genetically speaking, non-perfectness is recessive, and perfectness is dom- 
inant. 



2 There is one special case (m — 4 and 0(a) — abcb) where this is not true, due to 
the fact that B\ does not correspond to a ©2-word and because for this particular 8 the 
K-8-iuieLge of Sq equals S±. 
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Theorem 6. Let S and T be self- avoiding planefilling folding strings. If 
neither S nor T is perfect, then S * T is not perfect. If either S or T is 
perfect, then S *T is perfect. 

Proof: We use the characterization of perfectness based on the segment 
[1 + We write z s = f(9 s (a)), z T = f(9 T {a)), and z s * T = z s z T . 

Part 1: Suppose neither S nor T is perfect. Then K[9t(c)) + (1 + %)zt 
traverses [1 + %, i], and so part of the curve K[9s(9t(c))) + (1 + i)zs*r will go 
from the point (1 + i)zs to the point izs- But this part must be K[9s(c)] + 
(l+i)zs in the #s-loop, and since 9s is not perfect, this part traverses [1+i, i]. 
Conclusion: K[9s9t(c)]+izs*t traverses i], and thus S*T is not perfect. 
Part 2: a) First suppose S is perfect, T perfect or not. There is the special 
case that S = D, the 2-fold. Then 9s(a) = oh and 6*5(6) = cb, so 

9 s 9 T (a) = 9 s (ab . . . ) = abcb . . . 

starts with abc, and hence S * T is perfect. The following argument is illus- 
trated by Figure EJ If m s > 3, then either 9s(a) = abc . . . , and then [1 + i, i] 
is traversed by K[9s9t{o)}, or 9s{a) = aba . . . , in which case 1 + % is in the 
interior of the ^s^-loop, but also in the interior of the #s-loop (actually the 



(l + i)z s {l+i)z T 

, ^ , 




Figure 18: The 6>-loops of the non-perfect fold S = DUDDUUUDU, the perfect 
fold T = DUDD, and the perfect 50-fold S*T. 
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loop corresponding to the i^-^-image of the unit square, since the segments 
will not be traversed in the direct order). By perfectness of S, [1 + is 
traversed by either K[0g(ab)], which is the beginning part of K[9s9t(cl)], or 
by K[9s(d)] + izs which is the last part of K[9s9T(d)] + iz$*T- Therefore 
[1 + can not be traversed by K[9s9t{c)] + (1 + i)z s *T, and so S * T is 
perfect. 

b) This argument is illustrated by Figure [TBI If T is perfect, then (by 
Lemma [9]) part of one of the curves K[9 s (9 T (a))], K[9 s {9 T {b))} + zsz T or 
K[9s(0T(d))] + i z s z t will go from the point (1 + i)zs to the point i z$, and 
this part must be K[9s(c)] + (1 + i)zs in the #s-loop (or rather the K-9$- 
image of the unit square). Then if S is not perfect, this K[9s(c)] + (1 + i)zs 
traverses [1 +i, i]. So if T is perfect, and S not, then [1 + 2, i] is not traversed 
by K[9s9t(c)] + (1 + i)zs*r, and so S * T is perfect. □ 

The simplest example of a non-perfect fold, already given by Davis and 
Knuth Q3j), is the 4-fold S = DUU. Here K[9 n {a)\ fills a 45° - 45° - 90° 
triangle for all n, and the carrousel fills only half of the plane. 

Remarkably, the non-perfect 10-fold from Figure [TS] can be written as 
S = DUDDUUUDU = DUDD * U, and DUDD and U are both perfect 
folds. This does not contradict Theorem [61 since we always require S and T 
to start with D. The Q 2 of this perfect 50-fold is depicted in Figure [191 




Figure 19: The carrousel Q 2 of the 50-fold S * T of Figure [13 
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9 Boundaries 



Let w be a D 2 -word. The curve K[w] starts in 0, and ends in f(w). The 
smallest simple closed polygonal curve that contains K[w] and has vertices 
0, f{w) and further only mid grid points will be called the mid grid border 
of the curve X[u>], see Figure l20| Left. We have already encountered such 
borders: the octagon Q from Section [7] is the mid grid border of the curve 
K[6(a)} generated by any of the morphisms 9 that are constructed in the 
existence proof. 




Figure 20: Left: Mid grid border of the curve K[abcbcdcbcdadcdcb]. Right: 
the polyomino K [abcbcdcbcdadcdcb]. 

There is an alternative and attractive way to define this border, given 
by Knuth in the Addendum of p2] (see also the recent |13j): the mid grid 
border of K[w] is the topological border of the polyomino generated by w. 
In the context of Section [3] this polyomino is obtained via the diamond map 
K [-} mapping words to compact subsets of C, given by 

KM = B(l §), K.[b] = B(i, i), KM = B(-l i), K«[d\ = B(-l, |), 

where B(z, r) is the ball with center x and radius r in the Li-metric. Equation 
(JSJ) then defines the polyomino X<>M associated to any word w. 

Let S be a folding string with associated morphism 9. Gaussian integers 
that are visited twice by the curve K[9(a)] will be called internal points, and 
those visited once boundary points. Obviously, if S codes an m-fold, and we 
denote by i(S) and r(S) the number of internal and boundary points of S, 
then 

m = 2i(S)+r(S) - 1. 

Proposition 1. Let r(S) be the number of Gaussian integer boundary points 
of the curve associated to a fold S. Let S and T be the strings of two self- 
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avoiding folds. Suppose that S is planefilling and symmetric (S = S). Then 

r(S*T) = K(S)r(T). 

Proof: Let 9 s and 9t of length mj and my be the folding morphisms of 
S and T. It can be deduced from the symmetry that the mid grid border 
of K[9s{e)}, e = a,b,c,d consists of four congruent quarter-borders, each 
bordering say, ks boundary points (cf. Figure [21]) . Since S is planefilling, 
there are no other boundary points, except for the first and the last point 
of K[9s(e)}. Hence r(S) = 4k$ + 2. The curve K[6s0t{o)\ generated by 
S * T, is obtained by replacing the scaled segments of Ls(K[9t(cl)]) by 
translations of curves K[9s(e)}, e = a, b, c or d, where L$ is multiplication 
by z s = f(9 s (a)). 



Figure 21: The quarter-borders of the symmetric self-avoiding planefilling 25-fold 
9s (a) = abadadcbcdadcdadcbcdadaba in K[0s(c)] + 2zs in K[9s&t(o>)], where T is 
the 8- fold with 9t{o) = abadcdcd. In this example ks = 3. 

In this process we pass i(T) (scaled) internal points twice, where we lose 
eight times ks points, and four endpoints. Further we pass r(T) (scaled) 
boundary points once, where we lose (except at the begin and endpoint of 
Ls(K[9t(o)])) one endpoint and two times ks points. We thus obtain that 

r(S*T) = m T r(S) - i(T) (8k s + 4) - (r(T) -2) (2k s + 1) 

= (2i(T) + r(T) - l)r(S) - 2i(T)r(S) - (r(T) - 2)r(S)/2 

= r(S)r(T)/2. □ 

A slightly surprising consequence of Proposition [1] is that if S is any 
symmetric planefilling string with no internal points, then the fold S * S has 
(m — l) 2 /4 internal points. 
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10 Tiles 



By Theorem (2.4) of [I] the sequence of sets Lg n (K [6 n (abed)]) converges to 
a compact set 7s, which we call the tile generated by S. Now if S is a self- 
avoiding and planefilling m-fold, then K[6(abcd)] is maximally simple, and 
this implies that there exists a set of m Gaussian integers 2) such that 

K[9(abcd)] = (J (K[abcd] + o) . 

But then also for all n = 1,2,... 

K[6 n+1 (abcd)] = |J (lf[0 ft (a6cd)] 

and so for all n = 1,2,... 

L e L e n - l K[6 n+1 (abcd)] = [j (L g n K[6 n (abcd)] + o) . 

Taking the limit as n — >■ oo we obtain the self-similarity equation 

Le(T s )={J (Ts + D). (9) 

Applying Lq 1 on both sides we see that 7s is the unique attractor of a 
contracting iterated function systeiro 

A set T satisfying (jUJ) is called a self-affine tile . Such tiles have been 
extensively studied, mostly under the condition that Lg(G) C G and that 
the digit set D is a complete residue system for G/L$(G), where we write 
G for the Gausssian integers. In this case 2) is called a standard digit set. 
The authors of [S] and [TS] determine under the standard digit set condition 
the Hausdorff dimension of the boundary dTs of 7s , with a rather involved 
algorithm in both papers (see [7] for an even more involved algorithm that 
does not require the standard digit set assumption). 

For our folding tiles 7s it is obvious that they will not have a standard 
digit set when m is even (all digits will be even Gaussian integers). For the 

3 This connection between recurrent sets as in 0], and language directed IFS's, also 
called recurrent IFS's, or graph directed IFS's has been known for a long time, see [14] . 
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case m = 2 we will still obtain a tiling of the plane by a scaling, translation 
and rotation of the well known twin-dragon tile (it has been observed by 
Bandt ([2], p. 550) that this holds for any 2 element digit set {0,0}). 

For symmetric self-avoiding planefilling folds S we present a particularly 
simple formula for diniH(<97s), where dimn(-) denotes the Hausdorff dimen- 
sion of a set. (Actually [8] and p3] do apply in this case, as it is easy to show 
that for symmetric S the digit set will be standard.) 

Theorem 7. Let S be a self- avoiding planefilling symmetric m-fold with r 
boundary points. Let Ts be the tile generated by S with topological boundary 
dTs- Then 

, . 21og(r) -log(4) 

dim H (<97^) = — ^—r . 

log(m) 

Proof: The tile T s is a union of 4 sets: the limits of L- n (K[6 n (e)}) + /(e), 
where e = a,b, c, d. By Lemma El the boundaries of these sets all have the 
same Hausdorff dimension. So it follows with the usual counting and scaling 
argument (cf. [10]), that 

dimH(<97s) = lim 



iog(M 



where z is the endpoint of K[8(a)]. Since S is self-avoiding and planefilling, 
we have \z\ = \/m, and from Proposition [1] we obtain 

/l\n-l /r\ n 

r(5-)=( 5 ) r» = 2(-). 

This leads to the formula in the theorem. □ 

We give an example: for m = 5 we have the symmetric fold S = DDUU 
with 8(a) = abcba. Since there are no internal points, the boundary of the 
tile Ts has Hausdorff dimension 

21og(6)-log(4) = log(9) 
log(5) log(5) ' 



11 Other angles, more morphisms 

The strip of paper can be opened at other angles than 90°. In their paper 
([3]) Davis and Knuth observe that the fold S = DU, when opened at 60° 
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angles yields an interesting curve, which they coined the terdragon. They 
prove, using a ternary number system, that 6 infinite terdragons rotated 
over successive multiples of 60° pass exactly once through every edge of the 
triangular lattice. Actually, the theory developed in this paper can give a 
completely different proof of this property. We should now use a 3 letter 
alphabet, say {a, c, e}. The letters correspond to segments in the complex 
plane by the map /, now defined by 

/(o) = l, f(c)=u, f(e)=u 2 , withu:=e 2 ^ 3 . 

The anti-morphism r is as before, and the morphism a is now defined by 
a (a) = c, cr(c) = e, cr(e) = a. Folding morphisms again are those that com- 
mute with err, and, moreover, do not contain aa, cc or ee. The basic example 
is 9 with 9(a) = aca, which generates the terdragon. 

Defining the #-loop by the curve K[9(ace)}, and the anti-#-loop by the 
curve K[9(aec)}, then Theorem [1] remains true provided we replace 'the 9- 
loop' by 'the #-loop and the anti-^-loop'. The difference is caused by the 
different symmetries of the triangular lattice. 




Figure 22: The #-loop of a symmetric and an asymmetric planefilling 75-fold. 

Theorem [2] and [3] will remain true with the obvious definitions of self- 
avoiding and planefilling, but the proof of the analogue of the existence The- 
orem H] fails miserably. The analogue of Theorem [4] would be that there is a 
self-avoiding planefilling curve on the triangular lattice generated by a m-fold 
for every m > 1 that can be written as 

m = r 2 + s 2 — rs, r,sG Z + , 

since \z\ = \/r 2 + s 2 — rs is the norm of an Eisenstein integer z = r + su, 
and we should have \z\ 2 = m. This set consists of the so called Loschian 
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numbers larger than 1 

C := {3, 4, 7, 9, 12, 13, 16, 19, 21, . . . }. 

The reason that the proof fails is that not all paths in the directed triangular 
lattice correspond to allowed folding words (as, e.g., the word aa). In fact not 
only does the proof fail, but the result is also not true, since one can show that 
there is no self-avoiding planefilling m-fold on the triangular lattice when m 
is even. (The proof considers the union of the #-loop and the anti-#-loop: if 
we remove from this set the set K[8(a)} then one obtains a simple curve that 
(as an application of Pick's Theorem shows) has m — 1 Eisenstein integers 
in its interior. By symmetry of the midgrid boundaries of the the #-loop 
and the anti-#-loop, these integers have to be distributed evenly over the two 
loops, which is only possible if m — 1 is even.) 

Nevertheless, self-avoiding planefilling curves generated by folds opened 
at 60° abound. For small m these all have to be symmetric, but for large 
m they can be asymmetric (see e.g., Figure [221 right). An infinite family is 
given, for example, by all m-folds with m = 3k 2 , k = 1, 3, 5, 7, . . .. The folds 
here end in 2 k + ku, and trace a paralellogram-shaped midgrid border in a 
zig-zag way (see Figure 1221 left for the #-loop of this fold when k = 5.) Here 
the midgrid border is defined as in the 90° case, by a diamond map K , where 
now, for example, K [a] is the quadrangle with vertices 0, | + 1 v3i, ~ — | v3i, 
and 1. 

Curves on the hexagonal lattice are obtained when we unfold at 120° 
degree angles. We need a 6 letter alphabet, say {a, b, c, d, e, f}, which cor- 
responds to the 6 th roots of unity. Theorem [T] remains true, when we define 
the 6*-loop in the usual way. See Figure [23] for how the 3-fold S = DD gen- 
erates a self-avoiding curve on the hexagonal lattice. Actually the limiting 
fractal is an old friend: it is easily shown that the curves Lq n K[6 n (a)} are 
exponentially close in the Hausdorff metric to the usual approximants of the 
Sierpinski triangle, also known as the Sierpinski gasket. The notion of plane- 
filling is problematic in the hexagonal case, since the grid points all have 
degree 3. 

The concept of self-avoiding and planefilling curves generated by mor- 
phisms on the 4 letter alphabet can be generalized in an other direction. In 
fact, a slight generalization of Theorems [2] and [3] will hold for morphisms that 
preserve D2-words (see [6]). From that paper we have the following example 
of a non-uniform morphism that generates a self-avoiding and planefilling 
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Figure 23: Folding in three according to S = DD, and opening at 120° degree 
angles. Shown are the #-loop and the 5 th order curve. 

curve: 

6(a) = abad, 6(b) = ab, 6(c) = cb, 6(d) = cdcd. 
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